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ABSTRACT 


The  principle  of  reciprocity  is  employed  to  formulate  integral 
expressions  for  the  fields  of  antennas  confronted  by  obstacles.  The 
expressions  thus  formulated  are  particularly  useful  in  applying  a 
high-speed  computer  to  the  calculation  of  antenna -obstacle  inter¬ 
action.  The  integrals  encountered  require  surface  integrations  over 
the  antenna  aperture  or  on  a  surface  enclosing  the  obstacle. 

The  problem  of  computing  the  far-field  pattern  of  an  antenna- 
obstacle  system  is  treated  in  detail.  These  calculations  can  be 
completed  by  employing  the  antenna  "aperture  fields"  and  the  plane 
wave  fields  scattered  from  the  obstacle.  An  alternate  method  of 
calculation  requires  a  knowledge  of  the  fields  set  up  on  the  obstacle 
by  the  antenna. 

A  number  of  approximate  methods  for  finding  the  plane -wave 
scattered  fields  and  the  fields  on  the  obstacles  are  discussed.  Well- 
established  techniques  exist  for  obstacles  very  small  or  very  large 
in  terms  of  wavelength,  but  treatment  of  scatterers  in  the  inter¬ 
mediate  range,  the  resonance  region,  is  more  complex.  The  investi¬ 
gator  must  choose  the  method  to  fit  the  problem  in  order  to  proceed 
efficiently.  For  this  reason,  several  different  techniques  are 
discussed  which  can  be  programmed  for  high  speed  computation. 
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The  validity  of  the  expressions  formulated  is  checked  by  a 
comparison  of  experiment  and  calculation.  The  far-field  principal- 
plane  patterns  of  various  horns  in  the  presence  of  long  cylinders 
were  measured  and  calculated  by  means  of  these  expressions. 
Excellent  agreement  is  obtained  even  for  cylinder -to -horn  separations 
of  only  a  few  wavelengths.  The  calculations  employ  plane  wave  - 
infinite  cylinder  boundary  value  solutions  and  an  integration  over 
the  antenna  aperture. 

More  complex  obstacle  configurations,  a  cylinder  and  die¬ 
lectric  sheet  combination,  and  finally  a  pair  of  parallel  cylinders 
are  treated  successfully  by  approximate  methods. 
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CHAPTER  I 


INTRODUCTION 

As  the  frontiers  of  space  are  expanded  and  the  problems  en¬ 
countered  in  the  fields  of  communication,  guidance,  detection  and 
tracking  become  more  complex  due  to  increased  velocities  and  new 
environments,  the  problems  associated  with  antenna  system  uesign 
and  analysis  also  become  more  complex.  Antenna-obstacle  inter¬ 
action  is  such  a  problem. 

The  obstacle  problem  is  becoming  more  acute  with  respect  to 
the  performance  of  high-speed  vehicle  antenna  systems.  Infrared 
systems,  streamlined  radomes,  pitot  static  booms,  radome  supports, 
and  other  antennas,  these  all  act  as  obstacles  and  affect  antenna 
performance. 

On  the  ground,  the  increased  size  of  antennas  for  long-range 
communications  necessitates  the  use  of  very  large  radomes  and 
large  antenna  feeds.  These  and  their  associated  supports  act  as 
obstacles. 

This  paper  deals  with  the  calculation  of  the  effects  of  obstalces 
upon  antenna  system  patterns.  The  formulation  of  a  solution  is 
made  so  that  any  extensive  study  of  an  obstacle  or  class  of  obstacles 
can  be  programmed  on  a  high-speed  computer. 
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Finally,  an  experimental  and  theoretical  study  of  a  number  of 
antenna- obstacle  systems  is  made.  The  results  verify  the  theory 
and  indicate  that  calculations  are  practical. 


CHAPTER  II 


THE  THEORETICAL  FORMULATION  OF  THE 
ANTENNA- OBSTACLE  PROBLEM 

A.  THE  FORMULATION  OF  THE  BASIC  RELATION 

The  analysis  of  many  phases  of  an  electromagnetic  problem 
can  be  completed  by  invoking  the  principle  of  reciprocity.  In 
particular,  when  a  scatterer  is  placed  in  the  near  field  of  an  antenna 
system,  expressions  for  such  quantities  as  antenna  gain,  pattern, 
boresight  direction,  and  impedance^  can  be  formulated  by  use  of 
the  reciprocity  theorems.  The  expressions  thus  deduced  are  in  the 
form  of  surface  or  volume  integrals.  It  has  been  shown^  that  a 
variety  of  integral  expressions  are  available  depending  upon  the 
surface  chosen  and  the  test  source  used.  Thus  that  form  of  the 
solution  to  a  given  problem  should  be  chosen  which  yields  numerical 
results  most  readily.  As  the  problem  is  formulated  the  shape  and 
electrical  properties  of  the  scatterer  as  well  as  the  size  and  illumi¬ 
nation  of  the  antenna  will  be  seen  to  affect  this  choice. 

The  physical  and  mathematical  significance  of  the  "reaction" 
and  reciprocity  are  well  treated  in  the  literatu'-e;^*  *  however,  the 
theoretical  formulation  of  the  antenna-obstacle  problem  will  be 
developed  from  fundamentals  for  the  sake  of  clarity. 
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Consider  an  infinite,  homogeneous,  isotropic  medium  with  the 
parameters  p.  ^  and  in  which  are  located  an  antenna  (a  source  of 
electromagnetic  energy)  and  various  scatterers  (passive  bodies  with 
electrical  properties  differing  from  the  infinite  medium).  Antenna 
1  is  located  in  a  volume  Vj  enclosed  by  the  surface  Sj  ,  and  all  the 
scatterers  are  located  in  a  different  volume  Vg  enclosed  by  Sg. 
This  is  illustrated  in  Fig.  1. 


Mo,  «o 


Fig.  1.  Antenna  1  and  the  scatterer  in  free  space. 

The  electromagnetic  fields  of  any  antenna  in  Sj  can  be  related  by 
Maxwell's  equations  as  follows: 

(1)  V  X  Ej  =  -jwp  1  Hi  , 

and 

(Z)  ^7  X  Hi  = 


j  U)  €  1  El  +  Ji 
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In  these  equations  and  throughout  this  paper  the  fields  are  assumed 
to  vary  as  with  respect  to  time,  and  the  current  depends 
upon  the  source  located  in  .  Finally  and  Cj  can  be  replaced  by 
and  €q  respectively  everywhere  outside  of  Vj  and  Vg,  while 
inside  and  Vs  and  Cj  depend  upon  the  antenna  and  scatterers 
postulated.  Thus  for  all  regions  outside  and  Vg  the  equations 

(3)  V  xEi 
and 

(4)  V  X  Hi  =  El 

hold.  The  fields  Ei  and  Hi  are  members  of  the  set  of  all  fields 
whose  source  is  restricted  to  the  volume  Vi  and  which  exist  in  an 
infinite  medium  arbitrary  only  inside  the  volume  Vg  , 

Consider  next  a  second  antenna  (No.  2)  located  in  a  volume  V2 
surrounded  again  by  an  infinite,  homogeneous,  isotropic  medium 
(Ho*  ^q)  «xcept  the  volume  Vg  which  can  contain  any  arbitrary 
passive  bodies  as  illustrated  in  Fig.  2. 

The  electromagnetic  fields  of  antenna  2  can  be  related  by 
Maxwell's  equations  as  follows. 


(5) 


V  X  ^2  -  - j  U)  ^i2  _Hj  » 
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(8)  VxHj  =  juje^Ej 

The  fields  E2  and  H2  are  members  of  the  set  whose  source  is 
restricted  to  the  volume  V2  and  which  exist  in  an  infinite  medium 
arbitrary  only  in  a  volume  Vg  . 

By  using  relations  (3),  (4),  (7)  and  (8)  and  the  vector  identity 
V  •  (A  X  ^)  =  ^  •  V  X  A  -  A  •  7  X  it  is  found  that 


(9)  7  •  (El  XH2  -  E,  X  H,)  =  0 

everywhere  exterior  to  Vi ,  V2  and  Vg.  A  volume  integration  of 
Eq.  (9)  over  this  exterior  volume  yields 


(10) 


E2  X  Hi )  dV  *  0  . 


The  divergence  theorem  can  now  be  applied  to  Eq.  (10)  since  the 
fields  are  well-behaved  functions  in  and  on  the  boundary  of  V.  Thus 


(11) 


^7  •  (Ej  xHj  -E2  xHi)dV»^(Ei  xlj  -£2  *  Hi ) 


*  n  dSaO 


where  Sis  the  surface  enclosing  the  volume  V  and  n  is  the  unit 
normal  directed  out  of  V.  The  surface  S  is  composed  of  the  internal 
boundaries  Sj ,  S2  »  and  Sg  and  an  external  boundary  Sq.  Therefore 
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Eq.  (11)  can  be  rewritten  as: 

(12)  ^  (Ej  xHj  -E,  xHj)  •  ndS  +^(Ei  xHj  -E*  xHj)'ndS=0 

Si  +S2  +S  ®o 


The  external  surface  integral  is  discussed  in  Appendix  A  and  does 
not  contribute  to  Eq.  (12)  so  that 

(13)  J  (El  xHj  -  Ei  xHi )  ‘^dS  +  (Ej  xH*  -  E*  xHi )  •  ndS 

Si  Ss 


I 


(Ej  xHi  -  ^1  xH^  )  •  n  dS 


This  relation  can  now  be  used  to  solve  many  aspects  of  antenna- 
obstacle  problems. 


B.  THE  FIELDS  OF  AN  ANTENNA  NEAR  AN  OBSTACLE 

Equ^.tlon  (13)  can  be  used  to  compute  the  fields  of  an  antenna 
transmitting  in  the  presence  of  an  obstacle.  The  fields  Ei  and  Hi 
appearing  in  Eq.  (13)  are  the  fields  of  just  such  an  antemia 
(See  Fig.  1). 
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In  order  to  find^i  at  an  arbitrary  point  in  space  (XqiYq,  Zq)> 
it  is  observed  that  within  V2 


(14) 


V  *  (El  xHj  -  Ej  xHi )  » ju(ii2  *  liz 


"jw  (^2  ”  £(,)  El  •  ^2  “ji  ’ 

Now  let  antenna  2  be  an  infinitesimal  filament  of  electric  current 
located  at  Xq,  yg,  Zq  with  ^2  »  ^2  =  o»  ^o  ^2  •  Then 


(15)  7  .  (El  X  H2  -  E,  xHi)=  -  £2  ‘El 
in  V2  and 

(16)  y  ^  *  ( El  XH2  -E,  xHi)dV  =  -  J  Ji  •  El  dV 

V,  V, 


The  divergence  theorem  yields 


(17) 


I 


(^1  *ii2  "£*  *iii )  •  ndS  =  -  \  £2  •  El  dV  , 


-I 


and  since  £2  is  an  infinitesimal  current  filament 


(18) 


^  (El  XH2  -Ej  xHi)  •  ndS  =  -Ei  (x^,  y^,  Zq)  •  f  2  I2 
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A 

where  Ij  and  ij,  describe  the  strength  and  orientation  of  J2  re¬ 
spectively*  Finally,  combining  Eq.  (13)  and  (18)  yields 


(19) 


^  (El  xHj  -  E,  xH,)  •  ndS  +  (Ej  xHj  -  Ej  xHj 


ndS 


“  El  (xQ.yo,  Zo)‘^z  h 


The  unknown  field  Ei  (xo.yo»  *0)  computed  in  terms  of 

the  above  surface  integrals  when  an  infinitesimal  electric  current 

filament  is  used.  The  surface  integrations  are  equivalent  to  a 

solution  of  the  boundary  value  problem  in  that  the  necessary 

boundary  conditions  are  imposed  by  specifying  or  knowing  the  fields 

nxEi  and  nxHi  on  the  surfaces  Si  and  Ss*  The  use  of  a  current 

filament  to  formulate  Eq.  (19)  is  a  modification  of  the  work  of 

2 

Dr.  J.H. Richmond. 

In  solving  Eq.  (19)  only  the  component  of  ^1  parallel  to  the 
current  is  found;  of  course  finding  all  three  components  of  Ei 
simply  amounts  to  the  use  of  three  mutually  perpendicular  dipoles. 
Finally,  other  quantities  of  interest,  such  as  antenna  gain,  pattern, 
etc. ,  are  found  by  examining  the  form  of  JEi  (xp.  y©,  *0)* 
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C.  SOLVING  FOR  THE  FIELDS 

The  result  of  the  previous  section,  namely  Eq.  (19)»  permits 
a  solution  for  the  fields  of  an  arbitrary  antenna  in  the  presence  of 
an  obstacle  provided  that 

1 .  the  fields  Ej  and  are  known  on  S^  in  the 
presence  of  the  obstacle; 

2.  the  fields  Ej ,  Hj ,  iiz  known  on 

the  surface  Ss; 

3.  the  fields  Ej  and  jji  of  an  electric  current 
filament  at  (x^,  y^,  z^)  are  known  on  Sj  when 
radiating  into  an  environment  arbitrary  within 
the  volume  but  free  space  elsewhere; 

4.  the  dipole  moment,  \  Jj  dV  =  ^2  > 

known  when  the  current  radiates  into  this 
environment; 

5.  and  the  surface  integrals  over  Si  and  S^  can  be 
evaluated  in  closed  form  or  by  numerical  techniques. 

It  is  appropriate  to  discuss  each  of  these  provisions  in  more 
detail. 
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1 .  The  antenna  aperture  fields.  The  f J elds  and  Hj  on 

some  surface,  Sj ,  enclosing  the  antenna  under  investigation  are 
referred  to  as  the  ’’aperture  fields”  of  the  antenna*  In  order  to 
employ  Eq.  (19)  and  compute  the  far  field  of  the  antenna  these 
aperture  fields  must  be  known. 

During  free  space  operation  the  aperture  fields  are  assumed 
known  through  measurement,  design,  or  previous  theoretical 
investigation.  The  free-space  aperture  fields  thus  constitute  the 
starting  point  for  finding  the  aperture  fields  in  the  presence  of  the 
obstacle. 

The  fields  and  H|  on  certainly  depend  upon  the  environ¬ 
ment  exterior  to  and  thus  can  be  escpected  to  differ  markedly  not 

only  when  an  obstacle  is  introduced  but  when  it  is  reoriented  or 
relocated. 

The  free-space  field  of  the  antenna  is  incident  upon  the 
obstacle  and  is  scattered.  This  scattered  field  then  is  incident 
upon  the  antexina  structure  and  is  rescattered.  This  process  of 
reflection  and  rereflection  continues  ad  infinitum.  Thus,  the  fields 
and  Hj  of  an  antenna-obstacle  system  consist  of: 

1.  the  free-space  fields  of  the  antenna, 

2*  the  fields  scattered  from  the  obstacle. 
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3.  and  the  fields  scattered  from  the  antenna  structure. 
Richmond^  has  shown  that  the  second  of  these,  the  fields  scattered 
from  the  obstacle,  do  not  contribute  to  the  integration  on  S]  in 
Eq.  (19)  . 

In  Appendix  B  the  contribution  of  the  fields  scattered  from  the 
antenna  structure  to  this  integration  is  discussed.  ^  particular 
antenna- obstacle  geometry  is  treated  which  points  out  the  factors 
that  affect  the  magnitude  of  the  error  introduced  by  neglecting  this 
contribution.  Such  an  approximation  amounts  to  using  the  free- 
space  aperture  fields  in  evaluating  the  integral  on  S|  in  Eq.  (19). 

The  experimental  and  calculated  results  of  Chapter  IV  and 
the  discussion  in  Appendix  B  indicate  that  such  an  approximation 
is  very  accurate  except  when: 

1.  the  antenna  and  obstacle  are  very  close  together 
(separation  less  than  4\  for  horn  and  cylinder), 

Z.  the  obstacle  is  near  the  focal  point  of  an  antenna 
system  or  near  a  large  plane  reflector, 

3.  or  the  antenna  input  is  poorly  matched. 
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2.  The  fields  on  the  surface  S^.  In  order  to  evaluate  the 
surface  integral  over  Sg  in  Eq*  (19)  the  fields  Ej ,  E j  ,  and 
must  be  known  on  In  general  these  fields  are  unknown  on  that 
surface,  although  in  particular  cases  their  measurement  could  be 
justified.  Theoretical  evaluation  of  and  Hj  on  Sg  constitutes  a 
problem  that  is  as  difficult  as  solving  the  antenna- obstacle  problem 
itself,  since  it  is  a  solution  of  that  same  problem  on  a  designated 
surface. 

It  is  fortunate  that  there  exists  a  particular  choice  for  the 
environment  within  Vg  of  the  current  filament  J2  that  causes  the 
integration  over  Sg  to  vanish.  This  is  discussed  in  Appendix  C, 
and  it  will  suffice  here  to  say  that  by  employing  the  fields  of  £2 
the  presence  of  the  same  obstacle  as  for  the  antenna  under  investi¬ 
gation,  the  integral  vanishes  and  Eq.  (19)  becomes 

(20)  J  (El  XH2  -Ez  xHi)*ndS-Ei(Xo.yo,Zo)*  l2  I*  . 

Si 

A  second  possibility,  of  particular  value  when  the  free- 
space  field,  is  known,  permits  a  calculation  of  the  field 
scattered  from  the  obstacle,  ^f,  by  an  integration  over  the 
scattering  surface  Sg.  In  order  to  show  this  Eq.  (20)  is  rewritten  as 


ft 
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(21) 


I 


(El  xHj  -  Ej  xHj )  •  ^dS  =  (E?  +Ei)  •  llz  I* 


and  Eq.  (19)  can  be  written  as 


(22)  (EfxHj  -E,  xH®)‘^dS+j'  (EfxHj  -E,  xH?)*ndS 


o  A 

^  2  Iz  • 


Based  on  the  discussion  of  the  antenna  aperture  fields  in  the  last 
section,  it  is  true  for  many  antenna- obstacle  geometries  that 

(23)  J(Ei  xH,  -E,  xHi)*  ndS=  ^(E^xH*  xlj®)  dS  . 


Therefore  combining  Eq.  (21).  (22),  and  (23)  yields 


(24)  ^  (E°  X  Hj  -  E,  xl^)  •  'll  dS 


8  A 

B  *  h 


and  by  a  similar  reasoning  process 


(25) 


(El  X  H  °  -  E  °  X  Hi )  n  dS  ■  E?  •  If  ^  Ij 
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3.  The  fields  of  Jj  *  As  has  already  been  shown  it  is  con¬ 
venient  to  postulate  a  current  filament  located  at  (xoiyoi  Zq) 
order  to  calculate  the  electric  field  at  that  point  due  to  an  antenna- 
obstacle  system.  Furthermore,  application  of  Eq.  (20)  results  in 
a  solution  when  the  fields  of  presence  of  the  scatterer  are 

employed.  From  this  equation  it  is  seen  that  the  fields  E]  and  Hj 
must  be  known  on  a  surface  enclosing  the  antenna  of  interest  (Sx). 

In  order  to  evaluate  these  fields,  it  is  helpful  to  examine  the 
fields  of  a  current  £  identical  with  jJj  but  radiating  into  free  space. 
Such  a  current,  when  oriented  as  in  Fig.  3,  produces  the  fields 
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so  that  phase  fronts  propagate  in  a  direction  normal  to  JE  and  H  and 
the  field  vectors  are  perpendicular  and  related  by  the  free  space 
impedance,  Zq,  in  any  small  region  of  space. 

Thus  if  a  current  £2  ^  great  distance  from  the  obstacle,  a 

solution  to  the  plane  wave  scattering  problem  will  yield  E2  and  Hi  • 
Here  then  is  the  utility  of  Eq.  (20)  —  that  the  plane  wave  solutions 
to  scattering  problems  can  be  used  to  calculate  the  far  fields  of 
antenna> obstacle  systems. 

4,  The  dipole  moment  of  »  The  current  filament,  £2  »  was 
postulated  mathematically  and  thus  need  not  be  physically  realizable. 
It  is  possible  to  further  postulate  that  the  value  of  the  dipole  moment, 

dV  ■  ii  •  is  fixed  and  does  not  depend  upon  environment  or 

"v. 

location.  Such  a  current  filament  corresponds  to  an  infinitesimal 
electric  dipole  fed  by  an  infinitesimal  constant  current  generator. 

5,  Evaluating  the  surface  integrals.  The  formulation  of 
integral  expressions,  such  as  Eq.  (20),  (24)  and  (25)  for  the  field 
of  an  antenna- obstacle  system,  has  been  completed,  so  that 
numerical  computations  are  possible.  The  possibility  of  evaluating 
the  integrals  in  closed  form  exists,  but  in  many  practical  cases  the 
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aperture  fields  are  known  from  measurement  or  design,  so  that 
numerical  integration  is  more  appropriate.  The  use  of  high-speed 
computers  makes  numerical  techniques  even  more  feasible  for 
solving  specific  problems. 

Finally,  as  demonstrated  in  later  sections,  proper  choice  of 
integration  surfaces  simplifies  the  calculations  considerably. 


CHAPTER  III 


APPROXIMATE  SOLUTIONS  TO  THE 
ANTENNA- OBSTACLE  PROBLEM 

A.  GENERAL  REMARKS 

In  Chapter  II  a  number  of  distinct  integral  expressions  are 
derived  for  the  field  of  an  antenna  in  the  presence  of  an  obstacle, 
Eqs.  (19)>  (20),  (24)  and  (25).  These  are  by  no  means  the  only 
possible  expressions  for  the  fields,  but  they  are  adequate  for  the 
present  discussion.  There  are  a  few  general  remarks  to  be  made 
before  going  into  the  details  of  solution. 

First,  the  existence  of  high-speed  digital  computers  and 
their  tremendous  capabilities  allow  the  investigator  to  use  numeri¬ 
cal  techniques  on  a  large  class  of  problems  heretofore  impractical. 
The  integrals  encountered  in  this  antenna- obstacle  study  are 
members  of  this  class. 

Second,  the  choice  of  which  expression  to  employ  in  solving 
a  particular  problem  ultimately  is  determined  by  an  estimate  of 
the  calculation  time  of  each. 

Finally,  the  integral  expressions  deduced  to  facilitate  calcu¬ 
lation  of  antenna- obstacle  patterns  require  a  knowledge  of  either  the 
plane-wave  scattering  of  the  obstacle,  Eqs.  (20)  and  (24),  or  the 
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tangential  fields  induced  by  the  antenna  on  the  scattering  surface, 

Eq.  (25),  Solving  for  these  quantities  becomes  the  greatest 
problem  to  overcome. 

Exact  boundary  value  solutions  have  been  obtained  for  some 
obstacle  geometries.  The  best  known  of  these  are  the  sphere  and 
the  circular  cylinder  of  infinite  length.  Expressions  are  available 
for  the  plane- wave  scattering  of  such  obstacles,  but  in  many 
instances  their  usefulness  for  numerical  calculation  is  limited 
because  the  infinite  series  encountered  converge  slowly.  Thus  any 
useful  treatment  of  antenzui- obstacle  interaction  which  is  applicable 
to  analysis  of  an  arbitrary  obstacle  must  consider  approximate 
techniques  for  evaluating  these  fields. 

B.  OBSTACLES  SMALL  IN  TERMS  OF  WAVELENGTH 

Lord  Rayleigh^  treated  the  problem  of  electromagnetic  waves 
incident  upon  small  obstacles  before  1900,  The  success  of  his  work 
is  attested  by  the  fact  that  scattering  from  bodies  much  smaller  than 
a  wavelength  in  cross  section  is  today  called  "Rayleigh  scattering"  . 
His  method  consists  in  viewing  the  incident  field  as  uniform  in  the 
vicinity  of  the  obstacle  at  any  instant  of  time.  Thus  an  electric  and 
magnetic  charge  distribution  on  the  scatterer  can  be  computed  by 
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employing  the  techniques  of  static  electricity  and  magnetism*  The 
harmonic  nature  of  the  fields  then  causes  a  small  obstacle  to  scatter 
as  equivalent  electric  and  magnetic  dipoles.  The  strength  of  these 
dipoles  depends  upon  the  values  of  the  incident  fields  and  the  electri¬ 
cal  and  geometrical  properties  of  the  body.  Thus 
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where  symbolizes  an  electric  dipole  and  a  magnetic  dipole, 
H^(ro)  a-nd  ^(1*0)  incident  fields  at  the  center  (of  mass)  of 

the  body,  and  a  and  b  are  dyadics  depending  upon  the  geometric  and 
electrical  properties  of  the  small  obstacle. 

In  order  to  compute  the  fields  of  an  antenna  system  in  the 
presence  of  such  a  scatterer,  use  of  Eq.  (25)  appears  most 
practical.  This  means  that  the  antenna  near-field  in  free  space 
must  be  given  or  computed  at  the  obstacle  location.  Such  a  calcu¬ 
lation  is,  at  most,  time-consuming  when  the  aperture  fields  of  the 
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antenna  are  known;  the  large  the  aperture,  the  longer  the  calculation 
tinie.  The  scattered  field  of  the  equivalent  dipoles  Eq.  (31)  and  Eq. 

(32)  can  then  be  found  by  means  of  Eq.  (25).  Finally  the  fields  of 
the  antenna-obstacle  system  are  given  by  the  sum  of  the  free- space 
antenna  fields  and  these  scattered  fields.  Again,  the  free-space 
antenna  fields  at  most  involve  a  lengthy  calculation,  but  this  can  be 
readily  programmed. 

A  second  approxinoate  approach  to  interaction  between  an¬ 
tennas  and  very  small  obstacles  is  the  use  of  volume  polarization 
currents.^  This  technique  is  useful  because  the  scattered  fields 
from  a  dielectric  body  can  be  thought  of  as  origituitlng  from  a 
volume  distribution  of  electric  and  magnetic  current  radiating  into 
free  space.  These  currents  occupy  the  dielectric  volume  and  are 
given  by 

(33)  Jp=jw(e-€o)E 
and 

(34)  Kp  a  ju)  li  * 

where  ^  and  c  are  the  parameters  of  the  dielectric  obstacle.  In 
general  E  and  H,  the  fields  inside  the  dielectric,  are  unknown. 
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However,  for  small  dielectric  bodies  whose  electrical  properties 
(permittivity  e  and  permeability  |j. )  differ  only  slightly  from  free 
space,  these  fields  can  be  approximated  by  the  incident  field*  Thus 

(35)  Jp  =  j«o(€-€o) 
and 

(36)  Kp  =  ju)  (HL-p.^)  . 

The  size  of  such  an  obstacle  again  makes  Eq.  (25)  the  most  practical 
approach  to  calculating  antenna-obstacle  interaction* 

The  polarization  current  approach  is  limited  to  dielectric 
bodies  which  differ  only  slightly  from  free  space  but  does  offer  an 
advantage  in  that  it  is  applicable  to  a  greater  range  of  obstacle  sizes, 
even  to  those  normally  thought  of  as  lying  in  the  resonance  region* 

c.  L/  jE  obstacles 

The  most  common  large  obstacle  encountered  in  antenna 
system  design  and  analysis  is  the  radome.  In  general,  the  exact 
solution  for  a  plane  wave  incident  upon  an  infinite  homogeneous 
plane  sheet  can  be  used  in  conjunction  with  Eq.  (20)  to  calculate 
the  antenna  pattern*  Such  an  approach  is  approximate  in  that  any 
physical  radome  is  of  finite  size  and  in  most  cases  is  not  plane. 
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This  plane  wave-plane  sheet  approximation  gives  useful  results 
in  all  but  the  most  streamlined  radome  problems.  In  such  cases  ray 
tracing  can  be  attempted  in  the  region  of  greatest  curvature  or  the 
excitation  of  surface  waves  postulated  as  in  the  case  of  a  wedge 
radome.'  Richmond^  has  shown  good  agreement  between  measure¬ 
ment  and  calculation  for  an  antenna  aperture  only  partially  ob” 
structed  by  a  dielectric  sheet.  The  calculations  are  made  by  inte¬ 
grating  on  the  surface  of  the  sheet,  i.e.,  Eq.  (24)  is  used. 

The  analysis  of  inhomogeneous  radomes  as  well  as  antennas 
surrounded  by  plasmas  can  be  handled  by  means  of  plane  wave  theory 
and  Eq.  (20).  The  difficulty  here  stems  from  the  complexity  of  the 
problem  of  a  plane  wave  incident  on  an  inhomogeneous  plane  sheet. 
Recent  work  by  Richmond*®*  **  provides  efficient  numerical  tech¬ 
niques  for  calculating  transmission  and  reflection  coefficients  as 
well  as  the  fields  inside  such  inhomogeneous  layers. 

In  order  to  compute  antenna-obstacle  interaction  for  large 
conducting  obstacles  with  large  radii  of  curvature,  physical  optics 
is  used.  Physical  optics  consists  of  a  special  approximation  of 
the  current  Induced  on  the  surface  of  a  perfectly  conducting  obstacle. 
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In  the  shadow  region  the  induced  current  is  assumed  to  be  sero 
while  on  the  illuminated  portion  it  is  assumed  to  be 

(37)  J«2^xH^  , 

where  n  is  the  unit  normal  to  the  surface  and  is  the  incident 
magnetic  field.  Such  a  current  would  indeed  exist  on  an  infinite 
plane  conductor  with  unit  normal  n» 

In  order  to  make  use  of  this  physical  optics  approximation 
while  calculating  the  electric  field  of  an  antenna- obstacle  system, 
it  is  most  appropriate  to  employ  £q.  (25).  The  free-space  magnetic 
field  of  the  antenna  is  first  computed  on  a  surface  corresponding  to 
that  of  the  obstacle.  The  physical  optics  current,  J^,  set  up  by  this 
field  is  then  used  in  Eq.  (25)  to  compute  the  scattered  pattern. 

A  less  known  approximate  technique  which  stems  from  the 
equivalence  principle*  ^  yields  an  ejqpression  for  an  equivalent 
magnetic  current  on  the  illuminated  surface, 

(38)  K  =  2E^xn  , 

and  zero  on  the  shadow  portion.  This  approximation  is  subject  to 
the  same  limitations  on  obstacle  size  and  shape  as  is  the  physical 
optics  approximation.  A  comparison  of  the  two  techniques  is  made 
by  Harrington*  *  and  it  is  shown  that  the  approximations  are  of  the 
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same  order.  Thus,  in  general,  the  scattering  as  described  by  1^  is 
as  close  to  the  escact  value  as  that  given  by  the  physical  optics 
current  jJ  . 

D.  OBSTACLE  GEOMETRIES  IN  THE  RESONANCE  REGION 

A  class  of  obstacles  exists  for  which  neither  large  nor  small 
obstacle  approximations  are  at  all  successful.  The  surfaces  of 
these  bodies  are  characterized  by  regions  of  small  or  medium  radii 
of  curvature  and  their  circumference  is  on  the  order  of  one  to  ten 
wavelengths. 

The  approximate  techniques  discussed  in  the  last  two  sections 
are  so  successful  in  solving  the  large  and  small  obstacle  problems 
that  emphasis  in  recent  years  has  been  upon  developing  methods 
for  such  "resonant  region"  geometries.  The  researcher  thus  has 
a  choice  of  the  methods  developed  by  numerous  physicists,  engi¬ 
neers,  and  mathematicians  in  this  area.  Unfortunately,  a  solution 
by  any  of  these  methods  usually  requires  an  overwhelming  amount 
of  numerical  work.  For  this  reason  the  discussion  of  approximate 
techniques  in  the  resonance  region  is  limited  in  the  following  pages 
to  those  methods  most  readily  programmed  on  high-speed 


computers 
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1 .  Iteration  of  the  exact  solution.  In  principle  any  electro¬ 
magnetic  scattering  problem  can  be  solved  by  means  of  the  boundary 
conditions  and  solutions  to  the  wave  equation.  However,  the  wave 
equation  is  separable  in  only  a  few  coordinate  systems^**  ‘  *  and  the 
surface  of  the  scatterer  does  not  always  coincide  with  one  of  these 
coordinate  surfaces.  On  the  other  hand,  the  straightforward  nature 
of  the  boundary  value  solution  and  the  knowledge  that  an  exact 
solution  can  be  formulated  are  good  reasons  for  considering  this 
method. 

Employing  a  complete  set  of  solutions  to  the  wave  equation  and 
applying  the  boundary  conditions  yields  equations  for  the  unknown 
coefficients  of  this  set.  Whenever  the  scattering  surface  is  a  com¬ 
plete  coordinate  surface  in  which  the  wave  equation  is  separable, 
e.g,,  sphere,  infinite  cylinder,  infinite  plane,  etc.,  the  boundary 
conditions  and  the  orthogonality  of  this  set  of  solutions  serve  to 
determine  the  unknown  coefficients,  and  the  problem  is  solved. 
Scatterer 8  that  are  not  complete  coordinate  surfaces  yield  an  infinite 
system  of  equations  in  an  infinite  number  of  unknown  constants  (e.g., 
slotted  cylinder,  ‘  *  conducting  strips,*  ^  etc.). 

Such  systems  of  equations,  although  impractical  for  hand 
calculation,  can  be  solved  approximately  on  high-speed  computers. 
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Often  the  infinite  syetem  of  equations  can  be  approximated  successful- 
U 

ly  by  a  finite  set.  The  number  of  equations  and  unknowns  needed 
to  yield  a  good  approximate  solution  depends  upon  the  size  and  shape 
of  the  scatter er.  The  ability  of  high-speed  computers  to  handle 
matrix  problems  and  employ  iterative  techniques  is  the  reason  that 
this  approach  can  be  successful. 

This  method  can  be  applied  directly  to  the  antenna- obstacle 
problem  in  either  of  two  ways.  First,  the  plane  wave  scattering 
from  an  obstacle  is  calculated  by  choosing  and  as  plane  wave 
fields,  the  resultant  scattered  fields  (Ej  ,  H2  )  can  then  be  used  in 
the  antenna  aperture  integration  of  Eq.  (20).  This  method  offers  the 
advantage  that  many  different  antennas  can  be  analyzed  using  the  same 
scattering  information.  In  order  to  compute  a  complete  antenna- 
obstacle  pattern  this  method  becomes  complex,  for  the  scattering 
information  must  be  known  for  plane  waves  incident  from  many 
different  directions. 

An  alternative  method  is  to  compute  the  fields  scattered  from 
the  obstacle  using  the  antenna  near  field  as  the  incident  field.  The 
scattered  fields  of  the  antenna-obstacle  system  are  then  given  by 
Eq.  (25).  The  disadvantage  in  this  technique  is  that  the  whole 
procedure  must  be  repeated  for  each  new  antenna  encountered. 
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2.  Multipole  eaqiangton.  Recent  work  by  Kennaugh*^  discusses 
a  technique  applicable  to  scatterers  of  finite  surface  area.  The 
method  is  of  particular  value  in  computing  the  approximate  scattered 
fields  of  bodies  in  the  resonance  region. 

The  completeness  of  the  characteristic  solutions  to  the  wave 
equation  in  spherical  coordinates  nuikes  it  possible  to  write  the  fields 
scattered  by  any  body  in  the  form  of  an  infinite  series  of  spherical 
vector  wave  functions;  one  expansion  is  necessary  in  the  interior  of 
the  scatterer  while  another  is  necessary  in  the  external  region.  A 
solution  to  the  problem  involves  an  infinite  number  of  unknown 
constants. 

Rather  than  using  iterative  techniques,  Kennaugh  approximates 
the  scattered  field  by  a  finite  series  of  N  characteristic  solutions  to 
the  wave  equation  in  spherical  coordinates.  This  multipole  ex¬ 
pansion  is  then  made  to  fit  the  boundary  conditions  at  M  points  on  the 
surface.  The  fit  of  the  boundary  conditions  is  made  by  minimizing 
the  squared  error  of  the  tangential  electric  and  magnetic  fields  at 
the  M  points.  This  makes  it  possible  to  choose  any  number  of  points 
at  which  to  enforce  the  boundary  conditions. 
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The  flexibility  and  systematic  nature  of  this  technique  makes 
it  appealing,  especially  for  use  in  conjunction  with  a  high-speed 
computer.  It  is  limited  in  that  the  accuracy  of  the  approximate 
solution  is  unknown  for  most  problems.  However,  the  accuracy 
does  depend  upon  the  size  of  the  body,  the  number  of  points  chosen 
and  their  location. 

The  method  originally  proposed  by  Kennaugh  involved  mini¬ 
mizing  the  integral  of  the  squared  error  over  the  surface  of  the 
scatterer.  This  technique,  although  more  difficult  to  program, 
offers  the  advantage  that  for  some  scatterers  the  set  of  constants 
associated  with  the  multipole  esqAnsion  are  identical  with  the  exact 
values  and  are  independent  of  the  number  of  multipoles  employed. 

Consideration  of  how  to  compute  antenna- obstacle  interaction 
employing  the  results  of  this  technique  closely  parallels  the 
discussion  in  the  last  section.  However,  the  investigator's  control 
of  the  form  of  the  scattered  fields,  i.e.,  the  multipole  expansion, 
would  often  make  the  antenna  aperture  integration,  Eq.  (20),  the 
most  practical  method. 

3.  Variational  method.  The  variational  method  consists  in 
first  finding  a  stationary  eiq>ression  for  the  quantity  of  interest  in 
terms  of  another  quantity  related  to  the  problem.  Relating  this  to 
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the  problem  under  discus sioi.  means  that  a  stationary  e}q>ression 
for  the  scattered  field  must  be  found  in  terms  of  the  tangential  fields 
on  the  surface  of  the  scatterer* 

Such  a  stationary  expression  has  the  property  that  a  first 
order  perturbation  of  the  esqpressions  for  the  tangential  fields 
results  in  a  second  order  change  in  the  scattered  field.  The 
second  order  change  is  small  if  the  tangential  fields  are  close  to 
the  correct  value.  This  property  in  turn  means  that  the  investi¬ 
gation  can  approximate  the  unknown  tangential  fields  by  known 
fiinctlons  of  arbitrary  amplitude,  Cj^.  These  constants  can  be 
determined  from  the  condition  that 


(39) 
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for  all  n  when  |  E’  |  is  stationary. 

In  order  to  make  use  of  a  high-speed  computer  the  variational 
procedure  must  be  carried  out  systematically.  A  sero-order  calcu¬ 
lation  is  made  with  a  zero-order  trial  function  and  then,  by  adding 
higher-order  functions,  higher-order  calculations  are  made 

until  successive  calculations  yield  the  same  result.  The  spherical 
wave  functions  are  appropriate  as  the  various  order  trial  functions 
for  nearly  spherical  scatterers  while  the  cylindrical  wave  functions 
apply  in  the  case  of  a  long,  thin  body. 
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The  non-mathematically  inclined  often  experience  difficulty 
in  formulating  the  stationary  expression.  Rumsey^”  has  intro¬ 
duced  a  quantity  called  the  "reaction"  which,  when  understood, 
yields  stationary  expressions  more  readily  to  the  investigator. 

The  reaction  concept  gives  physical  significance  to  the  quantities 
involved  and  insight  as  to  why  the  approximation  is  stationary. 

4,  Approximate  current  method.  Many  of  the  resonant 
region  problems  involve  perfectly  conducting  scatterers.  In  such 
cases  the  scattered  fields  are  easily  found  when  the  current  flowing 
on  the  scatterer  is  known*  Thus  approximate  techniques  for  finding 
this  current  have  been  studied  in  detail  both  e^erimentally  and 
theoretically* 

V*A*  Fock^  ^  pioneered  in  the  area  of  formulating  an  ex¬ 
pression  for  the  current  induced  by  a  plane  wave  on  a  perfect  con¬ 
ductor*  His  work  concentrated  on  two  dimensional  problems  such 
as  a  plane  wave  Incident  normally  on  a  circular  or  elliptic  cylinder. 
He  formulated  an  expression  for  the  current  in  the  region  of  the 
shadow  boundary  by  replacing  this  portion  of  the  surface  with  a 
parabola  and  then  reasoned  from  physical  considerations  that  the 
wave  equation  could  be  approximated  by  a  'samplified  differential 
equation*  The  current  computed  for  TE  polarization  by  this  method 
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compares  favorably  with  the  theoretical  and  ejqperimental  values  in 
the  region  of  the  shadow  boundary  for  cylinder  diameters  greater 
than  a  wavelength.  The  method  has  been  extended  to  three  dimension¬ 
al  problems  and  to  TM  fields  by  Goodrich.  ^  ^ 

The  exact  expression  for  the  current  induced  by  a  plane  wave 
on  a  circular  cylinder  of  infinite  length  is  well  known  but  is  written 
as  an  Infinite  series.  For  cylinders  cf  radius  greater  than  a  wave¬ 
length  the  series  converges  slowly  so  that  extensive  work  has  been 
done  in  obtaining  asymptotic  formulas  for  the  current.  Wetzel^^ 
and  others  **  make  use  of  the  asymptotic  form  of  the  Hankel 
function  and  thus  obtain  expressions  identical  with  Fock's  for  the 
current  on  a  circular  cylinder  in  the  region  of  the  shadow  boundary. 
This  formula  also  gives  agreement  in  the  deep  shadow  region,  and 
by  a  higher  order  approximation  of  the  Hankel  function  Wetzel's 
work  handles  even  small  cylinders.  The  result  of  computing  the 
current  cn  elliptic  cylinders  in  terms  of  the  radius  of  curvature  at 
a  point  and  the  corresponding  current  on  a  circular  cylinder  of  that 
radius  has  met  with  success.  ** 

The  formula  of  Fock  and  the  asymptotic  formulas  of  other 
investigators  can  be  given  physical  significance  which  in  turn  aids 
in  the  investigation  of  other  problems.  The  currents  on  cylinders 
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and  apheras  can  be  considered  as  decaying  waves  which  travel 
around  the  conducting  surface  even  into  the  deep  shadow  region* 

26 

An  enlightening  discussion  of  this  subject  is  given  by  King  and  Wu. 

A  method  applicable  to  dielectric  bodies  as  well  as  to  perfect 
conductors  has  been  developed  by  Keller. This  method  hinges 
on  the  concept  of  diffracted  rays  and  the  formulation  of  a  geometric 
theory  of  diffraction.  The  diffracted  rays  are  shown  to  obey  a 
stationary  principle  (Fermat's  principle  for  reflected  and 
refracted  rays).  The  biggest  problem  associated  with  this  method 
is  its  adaptation  to  a  high-speed  digital  computer  program. 

Recent  work  at  The  Ohio  State  University  Antenna  Laboratory 
by  Peters  and  Thomas^*  indicates  progress  on  the  problem  of 
scattering  from  dielectric  bodies  in  the  resonance  region.  Their 
method  is  based  upon  geometric  optics  and  gives  good  results  for 
backscatter  from  dielectric  spheres. 


CHAPTER  IV 


ANTENNAS  CONFRONTED  BY  CYLINDERS 

A.  CYLINDRICAL  OBSTACLES 

One  of  the  most  common  obstructions  encountered  In  antexma 
design  is  the  cylinder.  The  feeds  of  large  reflectors  are  often 
supported  by  metal  cylinders*  refueling  and  Instrumentation  booms 
upset  the  performance  of  systems  located  in  the  nose  of  aircraft 
and  missiles,  and  recent  ^tudies^*  of  metal  space  frames  for 
radome  support  Indicate  that  metal  cylinders  may  become  an  Integral 
part  of  hypersonic  vehicle  antenna  systems.  These  uses  as  well  as 
the  mathematical  and  geometrical  simplicity  of  the  cylinder  make  it 
an  ideal  starting  point  for  antenna- obstacle  calculations. 

The  results  of  previous  sections  Indicate  that  the  fields  of  an 
ant*  .ma- obstacle  system  can  be  found  by 

1.  finding  the  fields  of  J^2  on  the  appropriate  surface, 

2.  deducing  the  aperture  fields  of  the  antenna,  and 

3.  performing  the  integration  Indicated  in  Eq. (2C)  or  (24). 

The  second  and  third  of  these  steps  vary  little  from  problem  to 
problem  and  in  most  cases  only  consume  a  great  deal  of  calculation 
time.  The  first  step  can  he  completed  either  in  an  exact  maxmer  or 
by  the  approximate  techniques  discussed  in  Chapter  III. 
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B.  BOUNDARY  VALUE  SOLUTION 

In  order  to  apply  £q.  (20)  the  fields  of  the  current  filament, 

J]  ,  must  be  known  when  it  radiates  in  the  presence  of  the  obstacle* 
These  fields  can  be  deduced  by  considering  a  plane  wave  of  appropri> 
ate  amplitude  and  phase  incident  upon  the  obstacle  when  the  obstacle 
and  current  are  very  far  apart,  as  discussed  in  Chapter  II.  Such  a 
solution  yields  the  far  field  of  an  antenna-obstacle  system  and  is  the 
goal  of  this  portion  of  the  paper. 

The  solution  for  a  plane  wave  incident  upon  a  conducting  right- 
circular  cylinder  is  well  known’*’  and  will  not  be  derived  here.  The 
coordinate  system  and  geometry  are  shown  in  Fig.  4.  For  an  inci¬ 
dent  TE  plane  wave  from  a  direction  4  ^  the  fields,  in  cylindri¬ 
cal  coordinates,  external  to  the  cylinder  are 

(40)  Ej  «  0 

Ep  (p  k/sin  ej  y  fW  '^"eo) 

^  ^  n^l 

.  jk  2  CO*  60 

+  (koP  »in  60)]  *in  n(*-So)  e-*  o  , 
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(42)  E^  =  Ei(-j)^  [4(koP 


naO 


+  (k^P  sin  eo)lco®  «^^o*  • 


(43)  *  ■  ?  y  *^“60)1 

I  \  Jk.z  cos  60 
cos  n(«-()>  )  e**  o  , 


naO 


(44) 


and 

(45) 


Hp  a  —  (-j  COS  60)  2  ®n^”tJn<^oP  *^*^60) 


n«0 


'(1)  1  Jk  z  cos  e, 

+  C„  H„  (k^p  sin  e©)!  cos  nio-o^)  e  o 


^  ^  E,  j  cos  eo 


^  2q  ('  kQ  sin  0Q 


go 

y  2j"  [J„(koP  sin  e^) 

nal 


+  Cj^Hj^(kQp  sin  ©q)]  Bin  n  (o-Oq)  e^^o®  ®o 


In  theje  expressions  Ej  is  the  magnitude  of  the  plane  wave,  the 
primes  indicate  derivatives  with  ^espect  to  argument  and  Oj^  ■  i 
for  n  a  0,  On  a  2  for  nsl,2,  3***.  If  the  cylinder  is  a  perfect 
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conductor  of  radius  a,  the  constants  are  given  by 


(46) 


-  Jn  (^o  eo) 

i<n  *  t  " '  '  '  ' 

^  ©o) 


If  the  cylinder  is  composed  of  homogeneous,  isotropic  dielectric 
material  (pi^ ,  )  of  radius  a,  and  the  plane  wave  is  incident 

normally,  e^,  ■  v  /  2,  the  constants  are  given  by®  * 


(47) 


JTr  J„(ko»Jn(  "JF,  Jl^r)  Jn^^o^) 

iffr  •^n(^o^j  Mr^r^^n  '^n^^o*’JT^r^)^  (^o*) 


C„  » 


Similarly  for  a  TM  incident  plane  wave, 
(48) 


Ha  *  0 


(49) 


H  .  £i  (1-L - ^  y  2rj”  [J„(koP  sineo) 

P  Zo  Vpko  sineo/  ° 


h\  Ik  Z  cos  Qn 

+  Dn^n  »(♦“♦(,)  «  ° 


J  *n-^"  [Jnt^oP  “"©o) 

naO 


+  Dn^n^^^oP  ©o)]  Ms'to)  * 
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(51) 


W 

Ea  ■  -E,  ^  [Jn^'^oP 


naO 


+  (koP  sin  So)]  cos  n(«-*o)  ® 


jko*  CO.  00 


(52) 


and 


00 


Ep  =  Ej  (-jco.  eo)^  ej”  si"  ©o^ 


naO 


+  ®o)l  “(ik-Oo)  • 


(53) 


Ej  j  CO.  eo 


go 

^  2nj“  [J„(koP  sin  60) 


p  ko  .in  eo  nal 
+  (^oP  sin  6  )]  sin  n{^-^^)  e^^o®  ®o 


If  the  cylinder  is  a  perfect  conductor  of  radius  a, 


(54) 


'n 


-Jn(^o»  si"  ©o) 

^n^^i'o*  "i"  ©o^ 


If  the  cylinder  is  composed  of  homogeneous,  isotropic  dielectric 
()jif  I  e^)  of  radius  a,  and  the  wave  is  incident  normally 

P^r  r^r^”  Pr 

Dj^a  j _ 

P7  ^  Jn(i^o®' J P'T® 


(55) 
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The  above  TE  and  TM  ejqpreaeione  constitute  a  complete  solution  of 
plane  wave  scattering  by  conducting  cylinders;  any  incident  plane 
wave  can  be  e}q>ressed  as  a  linear  combination  of  TE  and  TM  plane 
waves.  Moreover,  plane  wave  scattering  by  dielectric  coated  metal 
cylinders  as  well  as  by  concentric  dielectric  cylinders  is  given  by 
these  same  eicpressions  when  the  appropriate  form  of  and  is 
used’  ’  (found  by  satisfying  the  boundary  conditions  on  the  cylindri¬ 
cal  surfaces  and  restricted  to  normal  incidence). 

Equations  (40)  through  (55)  yield  the  fields  E2  and  Hz  neces¬ 
sary  for  evaluating  the  integral  of  Eq»  (20).  The  performance  of 
antenxias  confronted  by  cylinders  was  studied  by  calculating  and 
measuring  the  patterns  of  antenna- cylinder  combixuitions. 

1.  Calculation.  As  has  been  indicated  the  calculations  were 
made  employing  Eq.  (20)  which  requires  an  integration  over  the 
antenna  aperture.  Equation  (20)  was  used  rather  than  Eq.  (24)  or 
(25)  for  several  reasons.  First,  Eqs.  (24)  and  (25)  involve  inte¬ 
grating  over  the  surface  of  the  obstacle  and  for  infinite  or  very 
long  cylinders  this  was  not  desirable.  Second,  the  calculations  and 
measurements  were  carried  out  upon  a  number  of  different  antennas 
for  which  good  estimates  of  aperture  fields  were  available. 
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The  calculation  of  the  plane>wave  scattered  fields  over  the 
antenna  aperture  was  first  carried  out  by  approximating  the  appropri- 
ate  infinite  series, Eqs*  (40)  through  (55), by  a  finite  number  of  terms. 
For  the  cylinder  diameters  used,  d  <  2\ ,  ten  terms  of  the  series  were 
found  sufficient.  The  number  of  terms  needed  to  calculate  the  scatter¬ 
ed  fields  suggested  a  second  method  of  calculation;  the  asymptotic 
form  of  the  Hankel  functions  was  employed,  simplifying  the  series. 
This  asymptotic  approximation  ii 

(56)  (kopsin©^)-- 

for  k^p  sin  ©^  »  1,  and  k^p  sin  ©^  »  n.  The  successful  use  of 
Eq.  (56)  thus  depends  upon  the  distance  between  the  anteimas  and 
cylinder,  p  sin  ©q.  Use  of  Eq.  (56)  also  depends  upon  the  diameter 
of  the  cylinder  which  established  the  number,  n,  of  significant 
terms  necessary  to  evaluate  the  series. 

The  asymptotic  approximation  yields  the  bi- static  scattered 
pattern  of  the  cylinder.  This  has  been  programmed  at  The  Ohio 
State  University  Antenna  Laboratory  and  was  used  for  some  of  the 


-jkjjp  sin©o 
koP  ©o 


1? 


calculations 
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The  aperture  Integration  of  Eq.  (ZO)  was  carried  out  on  a  desk 
calculator.  Simpson's  rule  was  employed  and  an  estimate  of  the 
error  was  obtained  graphically  so  that  a  sufficiently  small  Increment 
would  be  chosen.  E3q>erimental  measurement  served  as  the  final 
check  of  the  accuracy  of  solution. 

2.  Eacperiment.  All  of  the  antenna-obstacle  measurements 
were  made  Indoors.  The  antenna- cylinder  system  was  mounted  on 
a  turntable  and  rotated  about  a  vertical  axis.  Antcxma  patterns 
were  recorded  on  a  logarithmic  recorder  which  was  driven 
synchronously  with  the  turntable.  A  crystal  detector(lN53)  opex^ 
ating  in  the  square-law  region  measured  the  received  signal  at  the 
antenna  terminals. 

The  overall  lexxgth  of  the  range  was  thirty  feet  and  the  frequen¬ 
cy  was  35  kmc  (X  m  0.338  inches).  These  figures  limited  the  size 
of  the  systems  measured  for  which  reliable  far-field  results  could 
be  observed.  Antenna  diameters  greater  than  twenty  wavelengths 
or  antenna- cylinder  spacings  exceeding  twenty  wavelengths  were 
not  considered  for  this  reason. 

Finally  finite  lengths  of  metal  and  dielectric  cylinders  were 
used  to  approximate  the  infinite  cylinders  for  which  the  calculations 
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were  made.  The  cylinders  were  long  enough  so  that  a  change  in 
length  caused  no  noticeable  change  In  the  measured  patterns. 

C.  RESULTS 

1.  The  small  horn  and  a  cylinder.  The  first  step  taken  In 
obtaining  a  comparison  between  theory  and  e]q>erlment  was  the 
measurement  of  the  free-space  E  and  H  plane  patterns  of  the  an¬ 
tenna.  The  assumed  aperture  fields  were  used  In  Eq.  (20)  to 
compute  these  free-space  patterns.  Measurement  and  calculation 
were  then  compared  to  test  the  accuracy  of  the  aperture  field 
assumption. 

A  one  and  one-half  wavelength  square  horn  fed  by  rectangular 
waveguide  excited  In  the  TEqj  mode  was  the  first  antenna  so  tested. 
Results  of  the  free-space  measurements  and  calculations  appear  In 
Figs.  5  and  6.  The  phase  of  the  aperture  fields  was  computed 
using  the  apex  of  the  horn  as  a  phase  center  while  the  amplitude 
distribution  was  assumed  to  be  that  of  the  TEq^  mode.  The  close 
agreement  with  experiment  indicated  that  these  assumptions  were 
appropriate. 

The  aperture  fields  and  the  plane  wave  fields  scattered  from 
cylinders  were  then  used  to  calculate  the  patterns  of  this  antenna 
confronted  by  a  perfectly  conducting  cylinder.  The  diameter  of  the 
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cylinder  was  chosen  to  be  the  same  as  the  width  of  the  antenna 
(d  =  1.  5\  ).  H  plane  measurements  were  made  with  the  axis  of  the 
cylinder  perpendicular  to  the  plane  of  measurement,  as  shown  in 
Fig.  7.  In  Figs.  7  through  10  the  H  plane  patterns  of  the  one  and 
one -half  wavelength  horn  are  shown  as  the  antenna -cylinder  spacing 
was  varied  from  four  and  one-half  wavelengths  to  twelve  wave¬ 
lengths.  The  close  agreement  illustrated  further  attests  to  the 
validity  of  the  aperture  field  assumptions. 

The  cylinder  was  reoriented  so  that  the  axis  was  perpendicular 
to  the  E  plane  as  shown  in  Fig.  11  and  an  E  plane  measurement  was 
made  employing  the  one  and  one -half  wavelength  horn.  Measurement 
and  calculation  are  compared  in  Fig.  11;  the  agreement  again  is 
excellent.  Finally  a  dielectric  cylinder  was  introduced  into  the  near 
field  of  the  horn  and  an  H  plane  measurement  was  made.  The  cylin¬ 
der  axis  was  perpendicular  to  the  H  plane,  the  diameter  of  the  cylin¬ 
der  was  one  and  one -half  wavelengths  and  the  dielectric  constant  two, 
€y  =  2.  0.  A  comparison  of  calculation  arid  experiment  appears  in 
Fig.  12. 

These  comparisions  of  measurement  and  calculation  served 
to  establish  the  validity  of  Eq.  (20).  The  hand  calculations  were 


lengthy  but  not  impractical. 
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Fig.  7.  The  H  plane  pattern  of  the  1.  5\  horn  confronted 
by  a  1.  5\  cylinder  at  a  distance  of  6  \  , 
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Fig.  8.  The  H  plane  pattern  of  the  1.  5\  horn  confronted 
by  a  1.  5X  cylinder  at  a  distance  of  4.  5X  . 
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Fig.  9.  The  H  plane  pattern  of  the  1.  5\  horn  confronted 
by  a  1.  5X  cylinder  at  a  distance  of  9X  . 
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1  Fig.  10.  The  H  plane  pattern  of  the  1,  5X  horn  confronted 

by  a  1.  5X  cylinder  at  a  distance  of  12X  . 
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Fig.  11.  The  E  plane  pattern  of  the  1.5\  horn  confronted 
by  a  1.  5\  cylinder  at  a  distance  of  6\  . 


Fig.  12.  The  H  plane  pattern  of  the  1,  5\  horn  confronted 
by  a  1.5X  dielectric  (e^  s  2,0)  cylinder  at  6\  . 
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2*  A  horn  confronted  by  a  cylinder  and  dielectric  iheet.  Further 
work  was  done  using  a  five  wavelength  square  horn  fed  by  the  TEqi 
rectangular  waveguide  mode.  A  phasjs-correctlng  lens  was  us>^d  in 
the  horn  aperture  designed  to  yield  a  constant  phase  distribution.  The 
first  step  in  making  antenna- obstacle  calculations  employing  this  an¬ 
tenna  again  was  a  comparison  of  measured  and  calculated  free  space 
performance.  Fig.  13.  A  constant-phase  TEqi  amplitude  distribution 
yielded  the  calculated  pattern,  which  closely  approximated  that 
measured. 

The  first  antenna-obstacle  system  treated  consisted  of  the 
perfectly  conducting  cylinder  (d  ■  1.  5K  )  oriented  perpendicular  to 
the  H  plane,  six  wavelengths  from  the  antenna.  The  measured  and 
calculated  H  plane  pattern  appears  in  Fig.  14.  The  error  present 
at  the  minima  of  this  pattern  is  attributed  to  phase  error  present  in 
the  free- space  aperture  fields  of  the  antenna. 

A  plane  dielectric  sheet  (d  s  0.97\  and  =  3.9)was  added  to 
this  antenna -cylinder  configuration  as  shown  in  Fig.  15  and  measure¬ 
ments  were  taken  as  the  sheet  was  moved  with  respect  to  both  the 
antenna  and  cylinder.  These  measurements  appear  in  Figs.  15  thru 
19.  Figure  15  corresponds  to  the  case  of  the  sheet  tangent  to  the  side 
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Fig.  14.  The  H  plane  pattern  of  the  5\  horn  confronted 
by  a  1.  5X  cylinder  at  a  distance  of  6\  . 

of  the  cylinder  farthest  from  the  antenna  while  Fig.  16  corresponds 
to  the  case  of  the  sheet  tangent  to  the  side  of  the  cylinder  nearest 
the  antenna. 

In  order  to  calculate  the  H-plane  pattern  of  such  a  system  an 
approximate  approach  was  used  to  deduce  the  plane-wave  scattered 
fields  of  the  cylinder- sheet  combination.  These  fields  were  then 
used  in  £q.  (20)  to  calculate  the  antenna  pattern. 


For  Field  Oecibele)  Fo,  Field  (Oecibele) 
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Fig.  15.  The  5\  horn  confronted  by  a  1.5X.  cylinder  at  6X. 

and  dielectric  sheet  (e,  =  3. 9>  d  s  0. 97\  )  at  6.7  5\  . 


Fig.  16.  The  5\  horn  confronted  by  a  1.  5\  cylinder  at  6\ 
and  dielectric  sheet  at  4. 30\  . 
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Fig,  19,  The  5\  horn  confronted  by  a  1.5\  cylinder  at 
6\  and  dielectric  sheet  at  14X  , 

A  plane  wave  of  unit  amplitude  incident  upon  the  dielectric 
sheet  is  transmitted  as  a  plane  wave  traveling  in  the  same  direction 
but  with  an  amplitude  equal  to  the  plane- wave  plane- sheet  trans¬ 
mission  coefficient,  If  such  a  plane  wave  should  then 

encounter  the  cylinder  it  would  be  scattered  according  to  Eqs.  (40) 
through  (55)  with  replaced  by  Tj_(e)E2  .  These  scattered  fields 
would  then  interact  with  the  dielectric  sheet  and  be  reflected. 
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These  last  reflections  were  neglected  In  calculating  the 
patterns  shown  in  Figs.  15  thru  19*  This  approximate  solution  is 
applicable  to  pattern  calculations  of  antennas  confronted  by  structur- 
ally  supported  radomes.  The  agreement  indicates  that  the  higher 
order  reflections  can  be  neglected  with  some  success,  particularly 
when  the  sheet  has  a  transmission  coefficient  of  nearly  unity. 

3.  A  large  ground  plane  antenna  and  a  cylinder.  The  calcu¬ 
lations  made  thus  far  have  been  based  upon  the  assumption  that  the 
contributing  portion  of  the  antenna  aperture  fields  are  unaffected  by 
the  obstacle  as  discussed  in  Chapter  II,  Section  C.  This  assumption 
is  not  always  valid. 

The  antenxia  chosen  to  illustrave  this  point  was  an  open  ended 
waveguide  in  a  large  ground  plane.  A  perfectly  conducting  cylinder 
(d  ■  1.  5X  )  six  wavelengths  from  the  aperture  was  oriented  with  its 
axis  perpendicular  to  the  H-plane.  The  measured  H-plane  pattern 
is  shown  in  Fig.  ZO. 

The  H-plane  pattern  was  computed  assuming  that  the  aperture 
fields  were  unaffected  by  the  obstacle,  Fig.  20.  These  calculations 
differed  from  the  measurements  appreciably  because  the  cylinder 
induced  large  currents  on  the  ground  plane. 
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Fig.  20.  Open  ended  waveguide  in  a  large  ground  plane 

confronted  by  1.5X.  cylinder  at  a  distance  of  6\  . 


In  order  to  take  into  account  this  type  of  interaction,  image 
theory  was  used.  The  plane  wave  scattered  fields  of  the  cylinder 
and  its  image  were  employed  in  Eq.  (20)  to  calculate  the  antenna- 
cylinder  pattern.  Much  better  agreement  was  obtained  (Fig.  20). 
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4,  Antenna  confronted  by  a  pair  of  parallel  cylinders.  The 
last  antenna- obstacle  system  treated  consisted  of  an  open  ended 
waveguide  and  a  pair  of  perfectly  conducting  parallel  cylinders.  The 
axes  of  the  cylinders  were  perpendicular  to  the  E  plane  of  the  an¬ 
tenna,  and  the  plane  of  these  cylinders  was  six  wavelengths  from 
the  waveguide  aperture.  Each  cylinder  was  of  the  same  diameter 
(d  =  1.  5X  ),  and  the  cylinder  axes  were  separated  by  slightly  over 
four  wavelengths  (s  s  4.3X)  as  shown  in  Fig.  21. 


Fig.  21.  Open  ended  waveguide-two  parallel  1.  5\ 
cylinders  separated  by  4.  3X  . 
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The  pattern  of  the  system  with  one  cylinder  removed  is  shown 
in  Fig.  22  along  with  calculations  made  by  means  of  Eq.  (20).  The 
pattern  of  the  system  with  the  other  cylinder  in  place  and  this  cylin¬ 
der  removed  would  be  a  mirror  image  of  Fig.  22  about  the  zero 
degrees  scan  axis. 


Fig.  22.  Open  ended  waveguide- one  cylinder  off  axis  20°, 

In  order  to  calculate  the  E-plane  pattern  of  the  system  with 
both  cylinders  present  by  means  of  Eq.  (20),  the  plane  wave  scatter¬ 
ing  by  the  pair  of  cylinders  must  be  known.  The  success  of  neglecting 
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cylinder- sheet  reflections  In  part  two  of  this  section  suggested  that 
a  superposition  of  scattering  by  each  cylinder  alone  would  be  a  good 
zero-order  approximation.  The  result  of  calculation  and  experiment 
Is  shown  in  Fig.  21.  This  approximation  Is  seen  to  describe  the 
general  shape  of  the  pattern  but  errors  as  great  as  two  decibels  are 
present.  Thus  a  higher  order  approximation  seems  necessary. 


% 


CHAPTER  V 


CONCLUSIONS 

In  calculating  the  effect  of  an  obstacle  upon  the  performance 
of  an  antenna  it  is  convenient  to  have  an  e:q>res8ion  for  the  parame¬ 
ter  of  interest  which  can  be  programmed  on  a  high-speed  computer. 
The  integral  expressions  for  the  electric  field  of  an  antenna-obstacle 
system  Eqs.  (20),  (24)  and  (25)  are  of  this  type.  These  e3q>ressions 
are  very  general  for  they  apply  to  any  antenna -obstacle  geometry 
and  can  be  modified  to  obtain  other  pertinent  parameters  such  as 
impedance,  gain  and  boresight  direction. 

The  possibility  exists  that  the  integrals  can  be  evaluated  in 
closed  from  but  in  the  vast  majority  of  instances  the  field  quantities 
in  the  integrands  will  necessitate  numerical  integration.  Often 
these  quantities  are  known  from  measurement,  design  or  as 
solutions  to  plane  wave  scattering  problems.  In  these  cases 
numerical  Integration  yields  a  useful  solution. 

A  choice  of  which  expression  to  employ  depends  on  the  geome¬ 
try  of  the  problem  and  what  is  known  about  the  antenna  and  obstacle. 
Most  often  the  antenna  aperture  fields  are  known  suggesting  use  of 
Eq.  (20).  Small  scatterers,  on  the  other  hand,  suggest  an  inte¬ 
gration  over  the  obstacle  surface  as  in  Eqs.  (24)  and  (25). 
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Use  of  any  of  these  expressions  depends  on  a  knowledge  of 
either  plane-wave  scattering  by  the  obstacle,  Eqs,  (20)  and  (24), 
or  scattering  of  the  antenna  near-fields  by  the  obstacle,  Eq.  (25). 
Since  these  quantities  are  known  exactly  for  only  a  few  obstacle 
geometries,  approximate  techniques  are  discussed.  Several  such 
techniques  are  discussed  in  Chapter  HI.  A  choice  of  which  to 
employ  depends  upon  the  obstacle  geometry.  This  is  particularly 
true  in  the  resonance  region. 

Use  of  Eq.  (20)  to  calculate  antenna- obstacle  systems  per¬ 
formance  was  highly  successful  in  Chapter  IV.  The  calculation  of 
the  patterns  of  an  antenna  confronted  by  a  cylinder  was  completed 
employing  approximate  antenna  aperture  fields  and  the  fields 
scattered  from  cylinders  by  a  plane  wave. 

It  was  found  that  the  asymptotic  form  of  the  plane-wave 
scattered  fields  could  be  used  in  nnaking  the  calculations,  for 
the  cylinder  sizes  used  and  the  antenna-to- cylinder  spacings 
treated. 

A  very  simple  approximation  yielded  good  results  for  plane 
sheet-cylinder  combinations  corresponding  to  structurally 
supported  radomes,  and  some  success  was  observed  in  treating 
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combinations  of  parallel  cylinders.  These  results  are  applicable 
to  the  design  and  analysis  of  supported  or  sectioned  radomes  as¬ 
sociated  with  ground-based  antennas  or  antennas  for  hypersonic 
vehicles. 


APPENDIX  A 


THE  EXTERNAL  BOUNDARY  INTEGRAL 


the  external  botindary  of  the  volume  V  dlscuased  in 
Chapter  II  and  shown  In  Fig*  2  can  be  taken  as  a  large  sphere 
enclosing  the  smaller  volumes  V^,  and  The  existence 
of  such  a  surface  necessarily  implies  that  the  sources  and  scatter- 
ers  are  finite  in  extent  and  are  finitely  distributed  within  V* 

Ballantine^^  discusses  the  integral  and  attempts  to  show  that 
it  equals  zero  on  an  Infinite  sphere* 

The  integral  can  be  shown  to  be  zero  as  a  consequence  of 
Sommerfeld's  radiation  condition  as  the  sphere  recedes  to  infinity. 

A  somewhat  lengthy  but  straightforward  proof  is  obtained  by 
writing  both  sets  of  fields,  H|  and  Ej  *  jii  •  expansions  of 
outgoing  waves  such  as  the  outgoing  spherical  wave  functions* 

Such  an  expansion  is  valid  on  and  exterior  to  Sq  when  it  contains 
all  sources  and  scatterers*  Maxwell's  equations  and  the  orthogo¬ 
nality  of  the  spherical  wave  functions  on  a  sphere  then  cause  the 
Integral  to  vanish.  The  assumption  of  only  outgoing  waves  depends 
upon  the  radiation  condition,  but  the  integration  can  be  performed 
on  a  finite  sphere  rather  than  a  sphere  at  infinity. 
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APPENDIX  B 


THE  APERTURE  FIELD  APPROXIMATION 


The  free- space  aperture  fields  of  an  antenna  are  perturbed 
when  an  obstacle  is  introduced  because  of  multiple  reflections 
between  the  obstacle  and  antenna  structure. 

Richmond  ’  has  shown  tha;,  only  outgoing  fields  contribute  to 
integration  on  a  surface  enclosing  the  antenna.  Therefore,  only  the 
reflections  from  the  antenna  need  be  considered  as  perturbations  of 
the  antenna  aperture  fields.  In  order  to  obtain  some  idea  of  the 
factors  that  effect  the  magnitude  of  these  reflected  fields  a  special 
case  is  treated. 

Consider  the  antenna- obstacle  system  shown  in  Fig.  23.  In 
order  to  simplify  the  analysis  the  antexuia  and  obstacle  are  chosen 
to  be  infinite  structures  in  the  z  direction  with  properties  inde¬ 
pendent  of  z,  furthermore  the  antenna  fields  are  assumed  polarized 
as  shown  in  Fig.  23  (H  transverse  to^).  These  restrictions  reduce 
the  problem  to  a  two  dimensional  one. 

The  center  of  the  antenna  aperture  coincides  with  the 
coordinate  origin  while  the  obstacle  is  located  at  (R^,  . 


b6 
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Structure 

Fig.  23.  The  two  dimensional  antenna- obstacle  system. 


Consider  the  case  for  which  is  so  large  that  the  obstacle 
lies  in  the  far  field  of  the  antenna.  The  free-space  antenna  field, 
in  the  far  field,  can  be  written  as 


(57) 


El  (R,  (ji) 


in  V 

z  £«|  {^9  (^ )  “  z 


g- jkoR 

~w~ 


PM 


When  this  field  is  incident  upon  the  obstacle  scattering  occurs. 
At  a  great  distance  from  the  obstacle  the  scattered  field,  , 
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can  be  written  as 

“jk  r!* 

(58)  E°  (S',  ,').■!  ?  °  A(,')  Ef  (R„, .  )  . 

and  $^are  cylindrical  coordinates  measured  from  the  center  of 
the  obstacle,  Fig.  23.  and  A(/ )  is  a  complex  reflection  coefficient 
dependent  upon  the  shape,  orientation  and  constitutive  parameters 
of  the  obstacle.  Evaluating  this  field  at  the  center  of  the  antenna 
aperture,  R  a  0,  yields 


(59) 


o  A  e“‘^^o®’o  o 

Eg  (Rq.  it  )  =  S  -  A(ir )  E?  (R^.  . 


^  (Rq,  TT )  encounters  the  antenna  structure  and  is  scattered.  This 


scattered  field,  at  the  obstacle,  can  be  written  as 

(60) 


1  A  1  A 

Ei(Ro.  4>o)  *  *  E,  (Rq,  »  z 


e~'^^o^o  o 


where  B(^)  is  a  complex  reflection  coefficient  dependent  upon  the 
antenna  structure,  input  impedance,  and  antenna-obstacle  orien¬ 
tation. 


This  process  of  reflection  and  rereflection  continues  ad 
infinitum,  and  the  field  incident  upon  the  antenna  can  be  found  by 
summing  up  the  infinite  number  of  reflections  from  the  obstacle. 
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The  total  scattered  field  incident  upon  the  antenna  is 


(61) 

so  that 

(62) 


and 

(63) 


Ei(Ro*’f)  *  E8(®-o*’^)  +Eg(Ro,1T)  ••• 


t  /\  e”"^^o^o  o  1 

E-(Ro»’f)*'* -  M-fr)  (e,  (R  ,  it) +e}  (R  ,  tt) 


+  Ef(Ro.ir)  •••) 


t  e"*^^o^o  /  6”^^o^o 

Eg  (Rq*  If )  ■  *  — ;== - A(it  )  (  —7= -  P(*o) 


rs; 

e-J''o'^o 


1^0 


e“J^o^^o 

+  -W  ■--P(»o)  A(Tr)  B(»p)  -  - 

4  Rft 


g-jkoRo  e~'^^o^^o 

+  -_ — P(eJ  a*(,t)B*(^J  - ^ -  ••• 


This  expression  is  a  geometric  series  which  can  be  summed  when 

e-j^o^Ro 


(64) 


A(tt)  B(0o) 


<  1  . 


which  is  true  for  large  Rq*  Thus,  the  total  field  incident  upon  the 
antenna  is 


E;  (Rq.  tt  )  *  z 


e-j^o^Rp 
A  Rp 


A(tt)  P(0 


1  -  )  R(»q)  e~j^^o^o 

R« 


(65) 
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Eq.  (65)  can  be  rewritten  as 


(66) 


e‘ (Ro.^)  -  a 


A(ir)  P(»J 
R  -  A(tt)  B(^  ) 


A  similar  summation  yields  the  magnetic  field  incident  upon  the 
antenna, 

t  (-Rxz)  A(7r  )?(♦.)  Vo 

(67)  Hg(Ro‘'^)  =  -  - - -  • 

^  R^  -  A(ir)  B(«^)  e"-^^Vo 

When  these  fields  are  incident  upon  the  antenna  they  are  reflected 
from  the  metal  and  dielectric  positions  of  the  structure  as  well  as 
the  antenna  terminals.  At  a  great  distance,  the  field  scattered 
from  the  antenna  can  be  written  as 

t  A  BU)  e‘'^^o^  t 

(68)  E,  (R,  ♦)  =  z  -111- -  E^(Ro«^) 

4r 

^  A(Tr)  P(«  )  B(^) 

*  z  - 1 -  , 

|r  (Rq-  A(it)  B(*o)  e-j2koRo  ) 

where  B{^)  depends  upon  the  antenna  structure  and  its  orientation 


with  respect  to  the  obstacle 
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The  total  field  incident  upon  the  obstacle  is 


(®-o*  *o^  “  —1  ^^•o’  (^o*  ^o^ 


Equation  (69)  can  be  written  out  in  detail  as 
(70)  El  (Rq.  »o)  =  ^  f  ltp._ 


A(Tr)  B(»o)  e~j^^o^o 
R,  -  A(it  )  B(*  )  e“«>21'o®-o 


E  /R  ,  )  .  ^ 

JEiJ  '*^0*  ^O*  "  *  P 

(R-  -  A(it)  BU  )  e  J^hoRo  ) 


At  great  distances  the  total  fields  scattered  from  the  obstacle  can 


be  written  as 


t  /  ,  A  A(/)  e'j^o^ 

(R»  e  )  ®  — -  El  (Rqi  ^q) 

IF 


jA(/)P(»o)  IFq 

"  (R^-  A(Tr)  B(«o)  ^ 


h‘ (R'  «')  = 


(R'x  z)  A(/) 

Zo  Ir^ 


El  (Rq*  ♦(,) 


I 
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Finally,  the  electric  field  of  the  antenna- obstacle  system  can  be 
written  as  the  sum  of  the  free-space  antenna  field,  the  total  field 
scattered  from  the  obstacle,  and  the  total  field  scattered  from  the 
antenna.  Thus 

(74)  El  (R,  s)  ■  Ef  (R,  e)  +  e{  (R,  s)  +  e‘  (r',  ) 


and 

(7  5) 


(R.S) 


A 
I  Z 


[Ef(R,^) 


+  El  (R»  + 


A(4.') 

- 


(V‘o>  (Ro-*o>!  • 


An  approximate  solution  of  this  antenna- obstacle  problem 
which  neglects  the  fields  scattered  from  the  antenna  structure 
(i.e.,  it  is  assumed  that  Ei  •  0)  can  be  written  as 

^  ^  A(0')e-j^o^' 

(76)  Ei(R,«)  *  z  [Ei(R,«)  + - - - E 

,R' 


1  (Ro*  * 


The  error,  W(R,  s),  due  to  this  approximation  is  given  by  the 
difference  of  Eq.(75)  and  (76).  Thus 


y\r  t  A(s') 

W(R,«)  -  z  [E\  (R,*)  +  - 


■jV 


eJ 


(Ro'«o)l 


(77) 
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and  Eq.  (68)  yields 


(78) 


W(R,*)  •'i 


A(Tr)  P(»J 

[R^-  A(ir)  B(*^)  ] 


/^(»)  e“'^^o^  ^  B(*  J  A(s')  ■•■  ^o) 

\  ^  "  K~l^ 

when  the  point  of  observation,  (R,  ^),  is  far  from  both  the  antenna 
and  obstacle. 

Examination  of  this  e3q>ression  for  the  error  introduced  by 
approximating  the  aperture  fields  by  the  free*  space  aperture  fields 
leads  to  the  following  conclusions* 

1.  The  error  depends  upon  the  magnitude  of  the  reflection 
coefficient  of  the  obstacle  in  a  direction  back  toward  the 
antenna,  A(w).  The  smaller  this  coefficient—  the  smaller 
the  error, 

2.  The  error  is  directly  proportional  to  the  relative 
illumination  of  the  obsUcle  by  the  antenna,  P(so). 

The  angular  position  of  the  obstacle  greatly  affects 
the  accuracy  of  the  aperture  field  assumption. 
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3*  The  error  depends  upon  the  separation  of  the  antenna 
and  obstacle,  Rg*  For  the  case  considered,  Rq  large, 
the  error  is  inversely  proportional  to  the  separation 
and  rapidly  becomes  very  small  in  comparison  with 
the  approximate  solution,  Eq.  (77). 

4.  The  error  depends  upon  the  reflection  coefficient  of 
the  antenna,  B(e)  and  B(eg)>  In  some  instances  the 
major  contribution  to  antenna  scattering  is  a  reflection 
at  the  antenna  terminals.  In  this  case 

(79)  B(«)  .  CP(«)  , 

where  C  is  a  constant  dependent  upon  the  amount  of 
mismatch  at  the  antenna  terminals.  In  other  cases 
(e«g.,  antenna  mounted  in  a  large  ground  plane)  B(s) 
and  P((ti)  are  vastly  different  functions. 

The  treatment  of  this  two  dimensional  problem  serves  to 
indicate  the  type  of  approximation  involved  in  employing  the  free- 
space  aperture  fields  of  an  antenna.  Each  antenna- obstacle  system 
encountered  should  be  analyzed  in  terms  of  the  above  factors. 

Such  an  analysis  will  often  indicate  when  the  approximation  is  not 
expected  to  be  accurate  and  why. 
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The  axialysis  can  be  extended  to  the  case  of  arbitrary  polari¬ 
zation  by  first  considering  TE  fields  and  the  analogous  reflection 
coefficients,  B^(^)  and  Arbitrary  polarizations  are  treated 

as  a  superposition  of  TE  and  TM  waves  so  that  the  reflection  coef¬ 
ficients  become  two  by  two  matrices.  Obstacles  in  the  near  field 
of  antennas  further  complicate  the  problem,  for  the  Incident  fields 
cannot  be  treated  as  plane  waves  and  the  use  of  reflection  coef¬ 
ficients  such  as  B(«)  and  A(s^)  is  no  longer  justifiable.  In  this  case 
the  use  of  an  infinite  scattering  matrix  that  relates  the  field 
scattered  from  a  body  to  the  incident  field  is  appropriate.  In  the 
two  dimensioxxal  case  the  scattering  matrices  should  be  referred 
to  cylindrical  coordinate  systems. 

Treatment  of  three  dimensional  antenna- obstacle  geometries 
is  complicated  by  the  addition  of  another  variable,  but  closely 
follows  the  two  dimensional  case.  In  the  far  field  two  by  two 
matrix  reflection  coefficients  can  be  defined  that  depend  upon  the 
shape,  orientation,  and  constitutive  parameters  of  the  scatterer 
while  in  the  near  field  infinite  scattering  matrices  should  be 


employed. 


APPENDIX  C 


ELIMINATION  OF  THE  INTEGRAL  OVER  S, 


In  formulating  Eq.  (19)  the  environment  of  antenna  2  was  left 
arbitrary  within  the  volume  Vg*  By  an  appropriate  choice  of  this 
environment  Eq.  (19)  is  simplified. 

Case  1.  The  scatterer  is  a  perfect  conductor. 

The  surface  integral  over  Sg  from  Eq.  (20)  can  be  rewritten 
as 


(80) 


I  ( El  X  -  Ea  X  Hi )  •  n  dS 
Ss 


■I  ^ 


X  El 


H, 


A 

-  n 


X  Ej 


dS 


The  surface  Sg  Is  now  chosen  to  coincide  with  the  perfectly  con¬ 
ducting  surface.  The  necessary  boundary  condition  on  E]  is  that 
n  x^i  ■  0  on  Sg  so  that 


(81) 


xHj  -Ej 


xHi) 


n  dS  a 


(n  X  Ej  ) 


Hi  dS  . 


Furthermore,  if  the  perfect  conductor  is  also  a  part  of  the  environ 
ment  of  antenna  2  then  nxE.  *0onS_  and 

(82)  f  (El  X  H,  -  E,  X  Hi)  •  n  dS  »  0 

-'S 


s 
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so  that  Eq.  ( 19  )  becomes 


(83) 


£ 


(£l  xHj  -  £,  xHi )  •  n  dS  =  (Xp,  y^,  z^)  •  <2  I2 


Case  2.  The  scatterer  is  isotropic  and  described  by  the  complex 
scalar  quantities  pi(x,  y,  z)  and  cj  ( x,  y,  z  ) . 

The  surface  integral  over  Sg  can  be  rewritten  as 


(84) 


(  “52  *  Si  )  * 


(EixHj  -E2  x^)dV 


and  by  using  a  vector  identity 


(85) 


(Ei  xHj  -  El  xHi  )•  n  dS  = 
®s 


VxEi  -  El*  VxHi  +E2  •  7x  Hi  - 


Hi  •  ^xE2  )  dV. 


Inside  Vg  ,  which  is  source  free, 


(86)  V  x  El  =  -j  w  Pi  (x,  y,  z  )  Hi  , 


(87) 


V  X  Hi  =  jw  €1  {x,y,  z  )  El 
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and 

(88) 

7  X  ^2  *  “jw  Ma  (*»  y*  *)  H2  1 

(89) 

VXH2  » jw  €2 <**yt*)  £2  . 

so  that 

(90) 

f  (El  XH2  -  E2  X  H,)  •  A  dS 

[H2  *  (-j‘*>P'l  Iil)-^l  *  (jw62E2  )  +E2  •  (ju>€,  El) 

-Hi  •  (-jw^l2  H2)1  dV 

or 

(91) 

{  (El  X  H2  -  El  *  Hi )  *  " 

=  \  -ju)  [ Hj  •  Hi  i  "  ^2  )  +  El  •  £2  2  ■  1 1 )  ]  dV  . 

‘V 
*  8 


By  choosing  ^2  =  M’l  €j  *  ej  at  all  points  within  the  integrand 
vanishes  identically  so  that 


(92) 


xHi) 


ft  dS  s  0 


and  Eq.  (19)  again  reduces  to  Eq.  (83). 


The  surface  integration  on  Sg  thus  can  be  eliminated  by 
choosing  the  same  environment  for  antenna  2  as  for  antenna  1 
outside  of  the  source  regions  Si  and  S2  • 
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